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(u  Is  the  radian  frequency  and  T Is  the  sample  Interval).  After  this 
change,  the  magnitude  and  phase  for  various  output-input  ratios  is  plotted 
as  a function  of  wT.  Linear  plots  and  Bode  plots  (especially  useful  for 
real-time  control  system  design)  are  presented.  Equations  are  given  for 
the  frequency  at  which  the  output  Is  0.707  of  the  input,  for  the  minimum 
of  the  output  to  input  ratio  and  for  the  overshoot  (amount  by  which  the 
output  exceeds  the  input). 


A synopsis  of  previously  available  statistical  properties  of  a-B-y  trackers 
is  also  included. 
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INTRODUCTION 


The  a-8-y  tracker  Is  a system  of  equations  which  use  estimates  of  position, 
velocity,  and  acceleration  to  predict  a new  (updated)  position,  a new 
velocity  and  a new  acceleration.  Then  the  error  between  the  measured 
position  and  the  predicted  position  is  used  to  improve  (correct)  the 
estimates.  The  parameter  y is  associated  with  acceleration  and  if  it  is 
set  to  zero  the  a-8-y  tracker  becomes  an  a-8  tracker.  In  a-8  trackers, 
acceleration  is  not  estimated  nor  used  in  the  predictor  equations. 

The  o-b-y  tracker  Is  especially  well  suited  for  computer  algorithms.  The 
o-b-y  algorithms  require  very  little  storage  and  a small  number  of  compu- 
tations and  are  well  suited  for  real-time  applications. 

The  Kalman  tracking  algorithms  require  much  more  storage  and  many  more 
computations.  However,  the  Kalman  algorithms  produce  estimates  with  mini- 
mum rms  error  estimates  given  non-statlonary  measurement  noise. 

The  standard  a-e-Y  tracker  minimizes  the  exponentially  weighted  squared 
error  between  the  data  and  parabolic  line  through  the  data. 

The  QD  based  a-B-y  tracker  minimizes  the  squared  error  between  the  data 
and  a parabolic  line  through  the  data— over  some  number  of  points  and 
subject  to  slope  and  intercept  conditions. 


THE  o-b-y  TRACKING  EQUATIONS 


The  a-e-Y  tracking  Equations  (1)  through  (6)  are: 
Corrector  (or  Estimator)  Equations: 


x = x + a(x  - x„) 
n np  n np 


v = v + 
n np 


+ (,n  - xnp> 


a = a 
n np 


(xn  “ V 


(1) 

(2) 


(3) 


Predictor  Equations: 


x 

v 

a 


np 

np 

np 


Vi  * ViT  * -r  v/ 


vn-l  * VlT 


n-1 


(4) 

(5) 

(6) 


where 

xn  = measured  quantity, 

x . v.  a„  = estimated  quantities, 
n n n 

xnp*  vnp*  anp  = predicted  quantities, 

xnl , v anl  = estimates  from  previous  step,  and 

T = step  size. 

The  Equations  (1)  through  (6)  are  in  a form  that  can  be  efficiently  pro- 
grammed on  a computer. 

Substituting  Equations  (5)  and  (6)  into  (2)  and  (3)  yields  another  form: 


xn  = 

xnp  * o(xn 

' xpp>  • 

(7) 

7n  = 

Vi  + Vi 

T * + (xn  - xnp>  • 

(8) 

*n  = 

Vi  * -Jr 

!xn  ' V ’ 

(9) 

xnp 

■ Vl  * T*n 

+ _L_  7 TZ 

-12  Vr 

(10) 

PERFORMANCE  MEASURES  Til 
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One  performance  measure  is  the  ratio  of  output  variance  to  input  variance 
with  noise  as  the  only  input.  These  ratios  are  defined  by 

nj'H.  r.  Simpson,  "Performance  Measures  and  Optimization  Condition  for  a 
Third  Order  Sampled  Data  Tracker,"  Transactions  on  Automatic  Control, 
April  1963,  pp,  182-183. 
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k /q)  , steady-state  variance  in  position  output 
xv  ' variance  in  raw  position  input 

k (n)  = steady-state  variance  in  velocity  output 
v'  ' ~ variance  in  raw  position  input 

k / q \ _ steady-state  variance  in-acceleration  output 
av  ' ~ variance  in  raw  position  input 


The  equations  for  these  variance  reduction  factors  are 

kx (0) = 28y(2o2  + 2$  - 3aB)  - ay2(4  - 2a  - 8)  , (11) 

ky(0)Y2T2  = 4g3Y  - 482y2  + 2y3(2  - a)  , (12) 

ka(0)Y2T4  - 48y3  . (13) 

where 

Y2  = yC«(28  + y)  -2y](4  - 2a  - 8)  . (14) 


2 

One  transient  performance  measure  D is  obtained  from  the  response  of  the 
tracker  to  a unit  step  inout  of  velocity.  Ideally,  the  position  output 
would  be  a ramp,  the  velocity  output  would  be  a constant,  and  the  accelera- 
tion output  would  be  zero  for  t > 0.  The  Dz  measures  are  the  sums  of  the 
squared  error  between  the  ideal  response  and  the  actual  response: 

D2Y2  = 2y(2  - a) (1  - a)2  , (15) 

D2Y2T2  = 2Y[a2(2  - a)  + 20(1  - a)]  + (Y3/8)(2  - a) 

-(y2/4)[2(4  - 3a) (2  - a)  + a8]  , (16) 

D2yV  = «y2(4  - 2a  - 8)  + (y3/2)(2  - a)  . (17) 

2 

The  final  performance  measure  A Is  obtained  from  the  response  of  the  tracker 
to  a unit  step  input  of  acceleration.  Ideally,  the  position  output  would 


3 


be  a parabola,  the  velocity  output  would  be  a ramp  and  the  acceleration 
output  would  be  a constant.  The  A2  measures  are  sums  of  the  squared  error 
between  the  ideal  response  and  the  actual  response: 

A2Y2  = a(4  - 2a  - e)(l  - a)2  , (18) 

A i 

A2Y2T2  = (a/4) (8  - 2a ) 2 ( 4 - 2a  - B)  + (y/8)(2  - a)(4  - 2a  - B)2  , (19) 

A2Y2T4  = aB2(4  - 2o  - b)  + (y/2)[2o2  - (2  - a)f}](4  - 2a  - &)  . (20) 


Finally,  Simpson  determined  that  all  these  performance  measures  are  minimized 
if 


2B  - a(a  + 3 + y/2)  = 0 . (21) 

Note  that  Equation  (21)  is  one  equation  in  three  unknowns  and  does  not 
completely  specify  the  design. 


"STANDARD"  a-6-y  TRACKER 


If  the  error  between  the  sampled  data  and  a parabolic  line  approximating 
the  data  is  squared,  weighted  exponentially,  and  minimized,  a,  B,  and  y 
are  found  [2]  In  terms  of  a parameter,  0.  Kahri'ias  calls  this  tracker 
"critically  damped." 

a - 1 - 03  (22) 


B = (1  - e)2(l  + 0)  (23) 

y = (1  - 0)3  (24) 


Morrison  [3]  calls  such  a tracker  "a  fading  memory  polynomial  filter  of 
degree  2."  Morrison  [4]  gives  the  variance  reduction  factors  in  terms  of  e. 

[2]  P.  J.  Kahri 1 as , Electronic  Scanning  Radar  Systems  Design  Handbook,  Artech 
House,  1976,  p.  215. 

[3]  Norman  Morrison,  Introduction  to  Sequential  Smoothing  and  Prediction, 
McGraw-Hill,  1969,  p.  516. 

[4]  Morrison,  p.  526. 
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k (0)  = ^ - e)  (19  + 24e  + 16e2  + 6e3  + e4) 
x (1  + e)5 

k U(0)T2  • ( 49  * 508  ■>  1382_) 

v (1  + e)5  2 

k (0)T4  = 6 ^ 
a 0 + e)5 

QD  TRACKER  [5] 


(25) 

(26) 

(27) 


i 


The  QD  tracker  minimizes  the  sum  of  the  squared  error  between  m data  points 
and  a polynomial  passing  through  the  data— subject  to  intercept  and  slope 
constraints.  For  a "second  order"  QD,  the  polynomial  is  a quadratic  and  [6] 


a = 


60m 


10m3  + 33m2  + 23m  - 6 


8 = 


2a 

m 


Y = 


2a 


(28) 

(29) 

(30) 


a-B  TRACKER 


The  system  of  Equations  (1)  through  (6)  can  be  converted  to  an  a-B  tracker 
by  setting  y = 0.  Then  all  of  Simpson's  results  reduce  to  those  given  by 
Benedict  and  Bordner  [7], 

If  the  error  between  the  sampled  data  and  a straight  line  approximating  the 
data  is  squared,  weighted  exponentially,  and  minimized,  a and  6 are  found 
(see  Reference  [2])  in  terms  of  e. 

£5]  W.  A.  McCool , "A  Matrix-Vector  Formulation  of  the  QD  Digital  Filter," 
Technical  Report,  WSMR,  New  Mexico,  December  1968. 

[6]  W.  L.  Shepherd,  Draft  Report,  Instrumentation  Directorate,  WSMR,  New 
Mexico. 

[7]  T.  R.  Benedict  and  G.  W.  Bordner,  "Synthesis  of  an  Optimal  Set  of  Radar 
Track  While  Scan  Smoothing  Equations,"  Transactions  on  Automatic  Control, 
July  1962,  pp.  27-32. 
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T 


a = 1 - e2 
b = (1  - e)2 

Morrison  (see  Reference  [3])  calls  such  a tracker  "a  fading  memory  polynomial 
of  degree  1."  He  gives  (see  Reference  [4])  the  variance  reduction  factors 
in  terms  of  0. 


k (0)  = — -~-9-3  (5  + 40  + 02)  (31) 

* d + e)3 

k (0)T2  = -211-<  «2) 

V (1  + 0)3 

Z TRANSFORMED  g-B-v  TRACKER 


Substituting  Equation  (10)  into  Equations  (7)  through  (9)  and  taking  the 
Z transform  yields 


(Z  + a - l)Xn  + (a  - l)TVn  + [(a  - l)/2]T2An  = aZXn 
BXn  + (Z  + B - l)TVn  + (8/2  - 1 )T2An  = 8ZXn  , 


Y*n  + yTVn  + {Z  + y/2  ' 1)T  An  = yUr\ 


(33) 

(34) 

(35) 


Writing  Equations  (33)  through  (35)  in  matrix  form  and  using  Cramer's  rule 
yields 


*n  _ Z[g(Z  - l)2  ♦ (b  * y/2)(Z  - I!  jud. 

U A * 


7n  _ Z[b (Z  - l)2  + y (Z  - 1)] 
X.  - Ta 


jn_s  MZ...-1)  2- 


T2a 


(36) 

(37) 

(38) 


6 


1 


Xng _ = (a  + B + y/2)(Z  - I)2  ± (B  + 3/2  v)(Z  - 1)  + y 
Xn 

where 

A - (Z  - I)3  + a(Z  - I)2  + (b  + y/2)Z(Z  - 1)  + YZ  . 

STABILITY 


(39) 


i 


(40) 

i 


For  an  a-B  tracker,  the  characteristic  equation  Is  [8] 
A * (1  - a)  + (a  + 6 - 2)Z  + Z2  . 


(41) 


Stability  requirements  are  [8] 
a > 0 , 

B > 0 

2a  + B < 4 


(43) 

(44) 


Stability  constraints  for  an  o-B-y  tracker  may  be  obtained  by  expanding 
Equation  (40)  and  applying  the  constraints  given  In  Jury  [9]: 


0 < a < 2 , (45) 
2a  + B < 4 , (46) 
aB  + ay/ 2 - y > 0 , (47) 
Y > 0 . (48) 


Note  that  Inequalities  (42),  (43),  and  (44)  Imply  Expression  (45),  and  that 
Expressions  (45),  (47),  and  (48)  Imply  that  B > 0 (Expression  (43)).  Finally, 
observe  that  setting  y ~ 0 reduces  the  o-b-y  expressions  to  the  a-B  results. 
For  the  o-b-y  tracker.  Relations  (45)  through  (48)  are  necessary  and  suffi- 
cient. Expression  (47)  may  be  replaced  by  the  more  convenient  sufficient 
(but  not  necessary)  condition  of  (49): 

aB  > y • (49) 

{!]  It.' L'.  Sonzales,  "Performance  Models  for  Range  and  Angle  Trackers  Ap- 
plicable to  Phased  Array  Instrumentation  Radar,"  Electronics  Division 
Technical  Memorandum  68-2,  Instrjnentatlon  Directorate,  WSMR,  New 
Mexico,  June  1968. 

[9]  E.  I.  Jury,  Theory  and  Application  of  the  Z-lransform,  John  Wiley, 

1964,  p.  93. 
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FREQUENCY  RESPONSE  CURVES 


To  find  the  frequency  response  for  any  one  of  the  transfer  functions  (36) 
through  (40),  replace  Z by  e+J“T  and  plot  the  magnitude  and  phase.  (Recall 
that  the  frequency  resporse  of  digital  filters  is  periodic  with  period  2w  in 
wT.  Only  a half  period  is  shown  in  the  plot.) 

A linear  plot  of  |Xn/Xn|  ard  phase  of  Xn/Xn  Is  shown  in  Figure  1.  The  same 
plot  on  a log- log  scale  Is  shown  In  Figure  2.  Note  that  the  curves  of 
Figure  2 are  Bode  plots  and  can  be  used  In  control  system  design.  A design 
example  for  using  Figure  2 follows. 

4 

DESIGN  EXAMPLE 

1.  Select  the  sample  time  for  a standard  (critically  damped)  a-e-y 
tracker  with  o = .1  If  the  tracker  must  pass  5 Hz. 

2.  Determine  the  attenuation  of  a 60  Hz  signal  for  this  sample  time. 

3.  Determine  the  frequency  of  the  first  maximum  response  above  the 
desired  pass  band. 

SOLUTION 

1.  Use  3 dB  down  as  the  edge  of  the  pass  band.  From  the  graph, 
x(.l)  is  3 dB  down  at  wT  * .135. 


T = * -00*3  second  . 

Thus,  the  signal  must  be  sampled  235  times  per  second  (or  faster)  if  5 Hz 
Is  to  be  attenuated  no  more  than  3 dB. 


2.  uT  * x .135  - 1.62  . 

from  the  graph,  x(.l)  at  «T  * 1.62  Is  about  23  dB  down. 

3.  The  first  undesired  maximum  occurs  at 

U - 2it  y ,Q4-  - 1452/second  , 
or 

f * 231  Hz  . 


WT4 


\ 


I 
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ao  hz 


If  the  23  dB  of  attenuation  at  60  Hz  Is  Insufficient  and  a critically 
damped  tracker  Is  desired,  then  a must  be  decreased.  Of  course,  decreasing 
a will  Increase  the  overshoot. 

If  the  phase  characteristic  of  Xn/Xn  Is  unsatisfactory,  then  the  phase 
characteristic  of  Xnn/Xn  should  be  examined.  Figure  3 reveals  that  the 
magnitude  curve  for  Xnp/Xn  Is  very  similar  to  the  magnitude  curve  for 
Xn/Xn,  but  the  phase  Is  monotonlcally  decreasing. 

Figure  4 exhibits  the  magnitude  and  phase  of  Xn/Xn  and  the  magnitude  and 
phase  of  200T2  ZTn/Xn  on  the  same  log-log  axes.  Note  that  the  magnitude 
characteristics  are  very  similar  for  wT  above  0.2. 

Figure  5 compares  the  QD,  ct-B,  and  a-B-y  trackers  on  the  same  log-log  axes. 
The  values  of  a,  e,  and  y used  In  these  cases  are  shown  In  Table  1. 


TABLE  1.  a-B-y  VALUES  FOR  FIGURE  5 


STANDARD  a-B-Y 

STANDARD  a-B 

go 

a 

0.3 

0.3 

0.3 

B 

0.0356 

0.0267 

0.0362 

Y 

0.0014 

0 

0.00219 

The  m 
m - 21 

required  to  produce  a ■ .3  Is 
values  and  Figure  6 compares 

16.6.  Table  2 shows  the  more  typical 
the  3 trackers  for  the  same  a. 

TABLE  2.  o-B-y  VALUES  FOR  FIGURE  6 

STANDARD  a-B-Y 

STANDARD  a-B 

go 

a 

0.246 

0.246 

0.246 

0 

0.0231 

0.0173 

0.0234 

Y 

0.00072 

0 

0.0011 
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ANALYTICAL  EXPRESSIONS 


The  plots  of  Figures  1 through  6 contain  all  the  frequency  domain  informa- 
tion. However,  analytic  expressions  for  the  minimum  and  maximum  of  |Xn/XQ| 
and  the  associated  frequencies,  and  the  frequency_at  which  |Xn/Xn|  = 0.707 
would  be  convenient  and  useful.  The  minimum  of  |Xn/Xn|  occurs  at  u»T  = ir  and 
is  given  by  Equation  (50): 

| -_H„  | = 2a  - 3 

1 X„  'min  4 - 2a  - 6 
n 

Note  that  Equation  (50)  is  independent  of  y. 

Unfortunately,  an  exact  solution  for  the  frequency  where  |Xn/Xn|  is  0.707 
would  require  the  symbolic  solution  of  a cubic  equation.  However,  setting 
| Xn/Xn | * 1 results  in  an  algebraically  tractable  expression.  The  exact 
solution  is  given  in  Equation  (51)  and  a good  approximation  is  given  in 
Equation  (52): 


uT  = 

arc  cos(l  - b)  , 

(51) 

o)T  » 

S2b  . 

(52) 

Now  observe  that  the  frequency  response  from  the  point  where  |Xn/Xn|  = 1 
to  the  minimum  is  nearly  a straight  line  on  the  log  curve.  Thus,  approxi- 
mating the  response  by  a straight  line  should  give  a reasonable  result. 
This  procedure  results  in 


(50)  | 


»T  - lO'-^lnla/fP  . (53) 

The  same  procedure  was  used  to  estimate  the  -6  dB  point  (-.15  in  the  exponent 
of  (53)  becomes  -.3)  with  the  results  shown  In  Table  3. 


TABLE  3.  APPROXIMATE  VERSUS  EXACT  UT  FOR  -6  DB 


a 

.1 

.3 

.488 

Approximate 

.194 

.65 

1.18 

Exact 

.196 

.65 

1.20 

i 


16 


An  exact  solution  for  wT  which  yields  a maximum  is  algebraically  intract- 
able. Examination  of  Figure  1 reveals  that  the  maximum  of  j Xn/Xn ) occurs 
near  the  midpoint  of  the  0 to  |Xn/Xn|  = 1 line  segment.  Thus,  u>T  is 
approximately  1/2  of  the  u>T  which  yields  |Xn/Xj  = 1. 

■ <7T  ■ _ (50 

Instead  of  calculating  the  maximum  of  | Xp/X^, { , it  is  more  convenient  and 
useful  to  find  the  overshoot.  The  oversnoot:  is  defined  overshoot  = 
maximum  - 1.  At  wT  , the  overshoot  is  approximately 

HlaX 

2 2 

overshoot  = -w-  + — — ^ — , (55) 

u 


where 


M = -|-  p2(l  - o)(n  + 6/3) 

n2  . P3q,4,ajjn2 


D = n2  + — 


16 

3 


8 


(56) 

(57) 

(58) 


wi  th 


Note  that  increasing  y decreases  n and  thus  increases  the  overshoot- 
provided  that  y < o6/2.  Table  4 compares  the  exact  overshoot  with  the 
approximation  of  Equation  (55). 


TABLE  4.  APPROXIMATE  AND  EXACT  OVERSHOOT 


a 

.1 

.3 

.488 

Exact 

Overshoot 

0.269 

0.232 

0.192 

Approximate 

Overshoot 

0.275 

0.261 

0.237 

The  approximate  and  exact  wTmax  are  compared  in  Table  5. 

TABLE  5.  APPROXIMATE  AND  EXACT  a)Tmax 


a 

.1 

.3 

.488 

Exact 

u»T 

max 

0.0402 

0.13325 

0.2440 

Approximate 

wT 

max 

0.0419 

0.13338 

0.2366 

SUMMARY  AND  CONCLUSIONS 


The  frequency  domain  response  of  an  a-B-y  tracker  provides  information 
which  is  useful  for  the  design  of  an  a-B-y  tracker.  The  Bode  plot  is 
very  helpful  if  the  a-B-y  tracker  is  to  be  incorporated  into  a control 
system. 

Equations  (50),  (53),  and  (55)  present  results  which  also  are  useful  in 
the  design  of  a-B-y  trackers. 
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